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A new approach to building models of generalized parton distributions (GPDs) is discussed that 
is based on the factorized DD (double distribution) Ansatz within the single-DD formalism. The 
latter was not used before, because reconstructing GPDs from the forward limit one should start 
in this case with a very singular function rather than with the usual parton density 

This results in a non-integrable singularity at /3 = exaggerated by the fact that /(/3)'s, on their 
own, have a singular Regge behavior for small /S. It is shown that the singularity is regulated 
within the GPD model of Szczepaniak et ah, in which the Regge behavior is implanted through 
a subtracted dispersion relation for the hadron-parton scattering amplitude. It is demonstrated 
that using proper softening of the quark-hadron vertices in the regions of large parton virtualities 
results in model GPDs H{x,^) that are finite and continuous at the "border point" x = ^. Using 
a simple input forward distribution, we illustrate implementation of the new approach for explicit 
construction of model GPDs. As a further development, a more general method of regulating the 
/3 = singularities is proposed that is based on the separation of the initial single DD f{/3, a) into the 
"plus" part and the D-term. It is demonstrated that the "DD+D" separation method 

allows to (re) derive GPD sum rules that relate the difference between the forward distribution 
f(x) = H(x,0) and the border function H(x,x) with the D-term function D(a). 



PACS numbers: 11.10.-z,12.38.-t,13.60.Fz 
I. INTRODUCTION 

The ongoing and future experimental studies of Gen- 
eralized Parton Distributions (GPDs) JHl] require theo- 
retical models for GPDs which satisfy several nontrivial 
requirements, such as polynomiality [5], positivity [6]-[8] 
hermiticity [1 , time reversal invariance [5], etc., following 
from the most general principles of quantum field the- 
ory. In particular, the polynomiality requirement, which 
states that the x'^ moment of a GPD H{x^ ^; t) is a poly- 
nomial in ^ of the order not higher than n + 1, is a con- 
sequence of the Lorentz invariance. The polynomiality 
condition is automatically satisfied when GPDs are con- 
structed from Double Distributions (DDs) ^[1110, 
(see also [l0]), thus the problem of constructing a model 
for a GPD converts into a problem of building a model 
for the relevant DD F(/3, a; t). 

Since a DD F(/3, a; t) has hybrid properties: it behaves 
like a usual parton distribution function (PDF) with re- 
spect to /3, as a meson distribution amplitude (DA) with 
respect to a, and as a form factor with respect to the 
invariant momentum transfer t, it was proposed [0 [TT] 
(in the simplified formal t = limit) to build a model 
DD F(/3, a) as a product of the usual PDF /(/3) and 
a profile function /i(/3, a) that has an a-shape of a me- 
son DA. This construction allows one to get an intuitive 
feeling about the shape of GPDs and their change with 
the change of the skewness parameter ^. It was noticed 
[12], however, that in the case of isosinglet GPDs, such 
a Factorized DD Ansatz (FDDA) does not produce the 
highest, (n + 1)^* power of ^ in the moment of i^(x, ^). 
To cure this problem, a "two-DD" parameterization was 
proposed [12], with the second DD G(/3, a) capable of 



generating, among others, the required ^"^^^ power. It 
was also proposed [12] to use a "DD plus D" parameter- 
ization in which the second DD G(/3, a) is reduced to a 
function D{a) of one variable, the D-term , that is solely 
responsible for the ^"^^^ contribution. The importance of 
the D-term and its physical interpretation was studied in 
further works (see Ref. [13 and references therein). 

Later, it was found out that it is still possible to write a 
"single-DD" parameterization [14] that incorporates just 
one function, but produces all the required powers up to 
^n+i ^]^ig representation also has a remarkable property 
that it allows, in principle, to invert the GPD/DD rela- 
tion, i.e., to obtain DD if GPD is known. So far, however, 
the single-DD representation was not used for building 
models for GPDs using the factorized DD Ansatz. The 
reason is that one should use much more singular func- 
tion f{l3)/j3 rather than just the usual PDF /(/3) for the 
GPD reconstruction from the forward limit. The combi- 
nation /(/3)//3, being an even function in the singlet case, 
has a non-integrable singularity at /3 = 0, even if f{j3) is 
finite at /3 = 0. Furthermore, the fact that PDFs /(/3) 
have a singular j3~^ Kegge behavior makes the problem 
even worse. 

In an independent development [TF, an attempt was 
made to implant the Regge behavior into a GPD model 
constructed in the spirit of the covariant parton model 
[16 , with the hadron-parton transition amplitude writ- 
ten in the dispersion relation representation capable of 
generating the desired Regge behavior through an ap- 
propriately chosen spectral density. To handle a > 0, the 
subtracted dispersion relation was used. The outcome 
was the claim [15] that the GPDs H{x^ ^) in this model 
have a singular {x — ^)~^ behavior in the vicinity of the 
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"border" point x = ^, which, if true, would ruin the ap- 
phcabihty of the perturbative QCD formahsm employing 
GPDs, since the latter works only when the GPDs are fi- 
nite and continuous across the border point x = ^. 

Our starting goal was to examine the model of Ref. [15] 
and to pinpoint the physical assumptions that resulted 
in the prediction of the singular {x — ^)~^ behavior (for 
an earlier analysis, see Ref. [17 ). As our analysis shows, 
the singularity follows from the use of the point-like ap- 
proximation for the hadron-parton vertices. For bound 
states, however, one expects that the hadron wave func- 
tion would generate an additional power-like (or even ex- 
ponential) suppression in the regions where the parton 
virtuality k'^ is large. We found, that if such a suppression 
is properly included in the model, the resulting GPDs are 
finite and continuous for x = ^. 

In our study, we also observed that the expression for 
GPDs derived from the model of Ref. [15 corresponds to 
a single-DD representation. Moreover, it has the struc- 
ture of a factorized DD Ansatz, but with the singularity 
at /3 = regularized by the subtraction made in the 
dispersion relation for the quark-hadron scattering am- 
plitude. Thus, the model of Ref. [15 (corrected for an 
appropriate softening of the hadron-parton vertices) gives 
a framework for building GPD models within the single- 
DD scheme. Using a simple, but rather realistic model 
for the input forward distribution (i.e., usual PDF), we 
illustrate, step by step, how to use this framework for the 
construction of GPDs. 

In particular, we found that the model produces a 
D-teim contribution, despite the fact that it uses only 
the forward distribution as an input. The formal rea- 
son is that the subtraction introduced in the dispersion 
relation differs from the subtraction that converts the 
original DD /(/3,a) into a (mathematical) "plus" distri- 
bution [/(/3,a)] + , which, by definition, cannot generate 
a D-teim. This observation raises the questions of a gen- 
eral nature about the separation of the D-teim from the 
initial DD /(/3, a) of the single-DD formalism. 

We found that the separation of /(/3, a) into the "plus" 
part [/(/3,a)]+ and the I^-term can be used to rederive 
the sum rule p!8H22] related to the dispersion relation 
for the real part of the DVCS amplitude [18 - 22 , and we 
also gave the derivation of another sum rule [19| proposed 
as the ^ ^ limit of that generic sum rule, and which 
relates the difference between the forward distribution 
f{x) = H{x^ 0) and the border function H{x, x) with the 
D-teim function D{a). 

The paper is organized as follows. To make it self- 
contained, we start, in Sect. II, with a short review 
of the basic facts about DDs and GPDs. In Sect. Ill, 
we describe the model [15^ with implanted Regge behav- 
ior, and give our derivation of expressions for GPDs and 
DDs that follow from this model. We stress the neces- 
sity of a profile function that eliminates the singularities 
for X = and present explicit results for models with 
two simplest non-fiat profiles. In Sect. IV, we perform 
a model-independent study of GPD sum rules, using the 




FIG. 1. Structure of the handbag diagram for deeply virtual 
Compton scattering. 

procedure of separating the initial DD into its "plus" 
part and the D-teim. We emphasize that H{x,0)/x and 
H{x^ x)/x^ due to their singular nature, should be treated 
as (mathematical) distributions rather than functions. 
Finally, we summarize the paper. 

II. PRELIMINARIES 
A. Double distributions 

Generalized parton distributions (GPDs) [T]-[4] natu- 
rally appear in the perturbative QCD description of the 
deeply virtual Compton scattering (DVCS) [23], [3] (for 
reviews see [5l [T3 | [24H27| ) , the process in which a highly 
virtual photon with momentum upon scattering on 
a hadron converts into a real photon with momentum 
q' = q -\- r. Basic features of GPD construction, in fact, 
are not specific to QCD, and may be illustrated on exam- 
ples of simpler theories [28 . In a toy scalar model (scalar 
quarks interacting with a scalar photon (/) through tptjjcj) 
vertex), the lowest-order (handbag) diagram (see Figjl]) 
may be written, in the coordinate representation, as 

C{q,P,r) = J p(^)e-^(^-)/2-^(^'-)/2 

x{P-r/2\ilj{z/2)ilj{-z/2)\P^r/2)d'^z , (1) 

where z is the separation between the "photon" vertices, 
and P = (p + p')/2 is the average of the initial p and the 
final momenta of the struck hadron, and V{z) is the 
quark propagator. 

The matrix element may depend on the coordinate 
difference z through invariants (Pz), (rz) and z'^ only. 
For large = —(7^, the higher terms of the z'^ expan- 
sion have suppression, thus the leading power term 
is generated from the matrix element taken at 2:^ = 0. 
The extraction of the = part of the matrix element 
may be performed in the standard way: through Taylor 
expansion in z followed by taking only the symmetric- 
traceless part (denoted by { }) 
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of the resulting local operators. For a scalar target, one 
may write 

{P + r/2\m{d^, ■ ■ ■ dMO)\P - r/2) 

oo r-n—l 



^=0 



Ml 



(2) 



In the momentum representation, the derivative d ^ con- 
verts into the average = {k^ + of the initial k 
and final k' quark momenta. After integration over /c, 
(kY should produce the P and r factors in the r.h.s. of 
the equation above. In this sense, one may treat (k)'^ 
as (/3P + ar/2)^ and define the double distribution (DD) 
[Il[3l[8l|9] 



(n-0 



dp da — Ar, 



(3) 



as a function whose P^~^a^ moments are proportional to 
the coefficients A^i . It can be shown [1] O [11] that the 
support region Vt is given by the rhombus \a\ + < 1. 
These definitions result in the "DD parameterization" 



{P-r/2\^{-z/2)^{z/2)\P^T/2) 
= [ F{(3,a)e-'^^^'^-'''^'^'^^^dpda^O{z^) . (4) 



of the matrix element. 



B. Generalized parton distributions 

Substituting the DD parameterization of the matrix 
element into the expression for the Compton amplitude, 
one obtains 



C{q,P,r)= [ F{l3,a)D 



/3P- 



d/3 da 



(5) 



where D[l] is the quark propagator in the momentum rep- 
resentation. Thus, the leading-twist term corresponds to 
a parton picture in which the initial quark carries mo- 
mentum /3P + (1 + a)r/2. Neglecting P^, (Pr) and r^, 
we get 

(^ + /3P+(l + a)r/2)2 = -Q2 + 2/3(Pg0 + (l + ce)N0 , 

i.e., /3 and a appear in the propagator in the combina- 
tion j3{Pq') + a{rq')/2 only. The latter may be written 
as x{Pq')^ with x = + where ^ = {rq')/2{Pq'). This 
redefinition leads to the parton picture in which the ini- 
tial quark carries momentum (x + ^)P. Introducing the 
generalized parton distribution (GPD) [H [2j [9] 



H{x,0 



F{P,a)5{x-(3-^a)d(3da , (6) 




FIG. 2. Parton picture in terms of DDs and in terms of GPDs. 



one can write the handbag contribution as 

T{q, P, r) = j'^ H{x, ^)D{q + {x + C)P) dx . 

One may try to define GPDs directly: 

{P + r/2\ip{-z/2)^{z/2)\P - r/2) 



e-i^(P^)H{x,0 dx + Oiz"^) 



(7) 



(8) 



However, an immediate question is what is the skewness 
^ in this definition? It cannot be treated as the ratio 
of {rz)/2 and (P^), since the ratio {rz)/2{Pz) cannot 
be the same for all points z. Hence, it is impossible to 
straightforwardly use such a definition in the expression 
([T]) involving a 4-dimensional integration over z. But, 
if one uses the DD parametrization and integrates over 
z^ then the scalar products (Pz) and (rz) convert into 
the scalar products {Pq') and {rq')^ respectively, since 
all other invariants, P^^r^^ {Pr) are neglected when they 
appear in the ratios with {Pq')^ {rq') or ^9 . In this 
sense, only the q' part of z is visible in the final result, 
and one may define GPDs by the formula ([8| in which z 
is substituted by a light-like vector n proportional to q' ^ 
say, by = q'^/2{Pq'). 

Still, the appearance of process-dependent quantities 
like (rq') and {Pq') in the definition of GPDs confronts 
the basic idea of the factorization approach that the par- 
ton distributions are process-independent functions. The 
standard "escape" is that {rq') / {Pq') in the GPD defi- 
nitions is substituted by an apparently "process-neutral" 
ratio r^ jP^ ^ supplied by information that P basically 
defines the "plus direction" and that some vector n de- 
fines the "minus direction" (for DVCS, n ^ q'). But this 
procedure creates a wrong impression that the definition 
of GPDs requires a reference to a particular frame. As 
shown above, one can define GPDs H{x^ ^) through for- 
mulas Q , ([6| , which do not refer to any particular frame 
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or process, and ^ is just some parameter. Of course, for 
each particular process, ^ should be adjusted to the kine- 
matics of the process, e.g., ^ = {rq') /2{Pq') for DVCS. 
Also, the parton interpretation of GPDs has the most 
natural form in the frame, where P for a massless hadron 
(and t = 0) defines the plus direction. 



C. D-term 

1. Scalar quarks 

Parameterizing the matrix element ([2|, one may wish 
to separate the Ann terms that are accompanied by ten- 
sors built from the momentum transfer vector r only 
(and, thus, invisible in the forward r = limit), and 
introduce the D-term [12] 



j D{a){a/2yda = Ar^ 



(9) 



as a function whose (a/2)^ moments give Ann- Within 
the DD-parameterization, the separation of the D-term 
can be made by simply using e~*^^^^^ = [e~*^^^^^ — 1] + 1. 
The D-term is then given by 

nl-\a\ 

D{a) = / F{P,a)dp , (10) 

and the DD-parameterization converts into a 
"DD plus D" parameterization 

{P-r/m-z/2)i;{z/2)\P + r/2) 

= [ [F{p, a)]+ e-'^iP-)-i<-i-r-)l-^ dp da 
Jn 



+ D{a)e-'°'^'''^/^da + 0{z^) , 



(11) 



where 



[F(/3, a)]+ = F(/3, a) - 5{(3) C F(7, a) dj . (12) 

is the DD with subtracted D-term. Mathematically, 
[F(/3,a)]+ is a "plus distribution" with respect to p. It 
satisfies the condition 



l-\a\ 



[F{f3,a)Ud(3 = , 



(13) 



guaranteeing that no D-term can be constructed from 



models with spin-1/2 quarks, it is more difficult to avoid 
an explicit introduction of extra functions producing a 
D-term. The basic reason [12] is that the matrix element 
of the bilocal operator in that case [30j should have two 
parts 



(P - r/2\^{-z/2)-f^^{z/2)\P + r/2)|twist-2 
= 2Pj{{Pz), {rz),z^) + r^g{{Pz), (rz), z^ 



(14) 



This suggests to introduce a parametrization with two 
DDs corresponding to / and g functions [12]. For the 
matrix element (14) multiplied by - which is exactly 



what one obtains doing the leading-twist factorization for 
the Compton amplitude [31] - this gives 

2^(P - r/2\i,{-z/2)j^i;{z/2)\P + r/2) 



-il3{Pz)-ia{rz)/2 



+ {rz)G{P,a) 



2{Pz)F{p,a) 
dp da +0{z^). 



(15) 



The separation into F- and G-parts in this case is not 
unique: expanding the exponential in powers of {Pz) 
and (rz), one may obtain the same {Pz)'^{rzy term both 
from the F-type and G-type parts. This leads to possi- 
bility of "gauge transformations" [32]: one can change 

F(/3, a) ^ F(/3, a) + 5x(/3, a) /da , (16) 
G(/3, a) ^ G(/3, a) - Sx(/3, a)/df3 , (17) 

using a gauge function xiP^ ^) that is odd in a. Still, the 
terms {Pz)^{rzy cannot be produced from the F-type 
contribution. The maximum of what can be done is to 
absorb all m / contributions into the F-type term. 
As a result, Eq. (15) is converted into a "DD plus D" 



parameterization |12 in which the term in the square 
brackets is substituted by the 



2{Pz)FD{i3,a)^{rz)6{(3)D{a) 



(18) 



combination, with D{a) given by the /3-integral of 
G(/3, a) and Fd{P^ a) related to the original DDs through 
the gauge transformation with 



XD(/3,a) = - / G(7,a)d7- / G(7,a) 



dj , 



(19) 



(cf. 



3. Spin-1/2 quarks: single-DD representation 



2. Spin-1/2 quarks: two-DD representation 

In the simple model with scalar quarks discussed 
above, one may just use the original DD F(/3, a) with- 
out splitting it into the "plus" part and the D-term. In 



In fact, since the Dirac index ja is symmetrized in 

the local twist-two operators ipij^ d ... d ^i^}"^ with 
the jii indices related to the derivatives, one may ex- 
pect that it also produces the factor PP^ + ar^/2. As 
shown by the authors of Ref. ^14j , this is precisely what 
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happens. In their construction, not only the expo- 
nential produces the 2;-dependence in the combination 
/3(Pz)+a(rz)/2, but also the pre-exponential terms come 
in the P{Pz) + a{rz)/2 combination, i.e., the result is a 
representation in which 



2iPz)F{p,a) + {rz)G{p,a) 
= [2l3{Pz)+a{rz)]f{l3,a) , 



(20) 



that corresponds to F{/3^a) = Pf{/3,a) and 
G(/3, a) = a). Thus, formally, one deals with 

just one DD f{P^a). In principle, though, this single 
function may be a sum of several components, e.g., 
5{a)f{l3)/l3 + 5{l3)D{a) /a (the result of the pioneering 
D-ieiia paper [12 for the pion DD in an effective chiral 
model corresponds to = 8{a)/\P\ — 8{P)/\a\). 

In the two-DD approach, GPDs are introduced 
through 

H{x,i)= I [F{p,a)^iG{p,a)\5{x- p -ia)dpda , 
Jq 

(21) 

which converts into 



(x,0= / {0^^a)f{0,a)6{x-0-^a) 

= x f{0^a)S{x — — ^a)d0da 
Jn 



(22) 



in the "single-DD" formulation. The D-term in the 
single-DD case is given by 



D{a) 



l-\a\ 



l + \a\ 



(23) 



D{a) 



a 



(24) 



and one may write f{0^a) as a sum 

f{P,a) = [f{(3,a)U+S{(3y- 
of its "plus" part 

a)U = f{0. a) - S{0) I /(7, a) dj (25) 



l-\a\ 



and D-term part 5{0)D{a)/a. 

Despite the fact that the "plus" part and the D-term 
are obtained from the same DD /(/3,a), they are inde- 
pendent in the sense that the "plus" part does not con- 
tribute in Eq. (23), and the i^-term contribution drops 
from Eq. (I25|). 



4. Getting GPDs from DDs 

The forward limit r = corresponds to ^ = 0, and 
GPD i^(x,^) converts into the usual parton distribution 
f{x). Using DDs, we may write 



f{x)= / F{x^a)da = x / f{x^a)da. (26) 

-'-l + l^cl J-l+\x\ 



Thus, the forward distributions f{x) are obtained by in- 
tegrating DDs over vertical lines = x in the (/3,a) 
plane. For nonzero ^, GPDs are obtained from DDs 
through integrating them along the lines = x — 
having 1/^ slope, i.e. the family of i7(x,<f) functions for 
different values of ^ is obtained by "scanning" the same 
DD at different angles. 



a 





























/3 



FIG. 3. Support region for double distributions and lines 
producing /(x),if(x,^) (for x > ^ and x < ^), if(^,^) and 

i^(-e,0- 

For X > ^ > 0, the integration lines lie completely 
inside the right half of the rhombus. The line producing 
GPD at the "border" point x = ^ starts at its upper 
corner, while the lines corresponding to \x\ < ^ cross 
the line = 0. Thus, one deals with the "outer" regions 
X > ^ and X < (in this case, the whole line is in the left 
half of the rhombus) and the central region — ^ < ^ < ^, 
when the integration lines in the (/3, a) plane lie in both 
halves of the rhombus and intersect the = line. 

In GPD variables (x, ^), the momentum fraction x — ^ 
carried by the final quark is positive for the right outer 
region, and negative for the central region, i.e., in the lat- 
ter case it should be interpreted as an outgoing antiquark 
rather than incoming quark [3j, i.e. GPD in the central 
region describes emission of a quark- antiquark pair with 
total plus momentum r+ shared in fractions (1 + x/^)/2 
and (1 — x/^)/2, like in a meson distribution amplitude. 

From this physical interpretation, one may expect that 
the behavior of a GPD H{x, ^) in the central region is 
unrelated to that in the outer region. But, since the 
GPD in both regions is obtained from the same DD, one 
may expect, to the contrary, that the set of GPDs for 
all "outer" x's and all ^'s contains the same information 
as the set of GPDs for all central x's and all ^'s. This 
"holographic" picture (cf. [20, 34 ) may be violated by 
terms contributing to GPDs in the central region and 
not contributing to GPDs in the outer regions: the terms 
with support on the = line, i.e., those proportional to 
S{0) (and, in principle, its derivatives), in particular, the 
D-teim. For this reason, the usual approach is to build 
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separate models for the i^-term and for the remaining 
part of DD. 



III. GPD MODEL WITH IMPLANTED REGGE 
BEHAVIOR 



Formulation 



5. Factorized DD Ansatz 



The reduction formula ( 26 ) suggests a model 



m 



h{P^ a) 



(27) 



where is the forward distribution, while h{f3,a) de- 
termines DD profile in the a direction and satisfies the 
normalization condition 



1-1/31 



1+1/31 



h{P^ a) da = 1 



(28) 



Since the plus component of the momentum transfer r 
is shared between the quarks in fractions (1 + a)/2 and 
(1 — a)/2, like in a meson distribution amplitude, it was 
proposed [8l[TT] to model the shape of the profile function 

by 



/ijv(/3,a) 



(1- 



2iAr 



\2N-\-l 



(29) 



with N being a parameter governing the width of the 
profile. 

Such a factorized DD Ansatz (FDDA) was originally 
[HI [11] applied to an analog of the a) function of 
the two-DD formalism, which corresponds to a model 
F{p,a) = f{j3)h{l3,a) and G(/3, a) = 0. Later, it was 
corrected by addition of the D-term |12j , which formally 
corresponds to the "gauge" (19) in which G(/3,a) 
Gd{P,ol) = S{(3)D{a), and F(P,a) Foi^.a). Note 
that \i F = j3f and G = af, the model FD{P^a) = 
h{l3^a) does not coincide with the model f{P^a) = 
f{P)h{P^a)/P^ since the gauge function xd{P^<^) (see 
Eq. ( |16| )) is nontrivial. 

Thus, there is a question whether the FDDA should 
be applied to F]j{P^a) (as it was done so far) or to 
the DD f{p^ a) of the single-DD formulation. It should 
be confessed that no enthusiasm has been observed to 



use FDDA in the form of the single-DD formula (27). 



This observation has a simple explanation: the function 
/(/3)//3 is not integrable for /3 = 0, even if f{/3) is finite 
for /3 = 0. The reason is that the DVCS amplitude con- 
tains singlet GPDs, which are odd functions of p. Hence, 
f{/3)/l3 should be an even function, and the principal 
value prescription does not work. Moreover, for small /3 
one would expect that the forward distribution /(/3) has 
a singular /(/3) ~ l//3^ Regge behavior, which makes the 
problem even worse. 



The assumptions used in the factorized DD Ansatz are 
based on the experience with calculating DDs for trian- 
gle diagrams [28^ and form factors in the light-front for- 
malism models with power-law dependence of the wave 
function on transverse momentum (see also [36 ). 

The simplest triangle diagram (see FigQ in the scalar 
model corresponding to Eq. ([2| may be used as an exam- 
ple of a model for GPD 



I 



d^kd{x-{kn)/{Pn)) 



{m\ - k\){ml - kl)(m\ - {P - k^) 



(30) 



Though the ^-dependence is not immediately visible 
here, it appears after integration over k through the 
{rn)/2{Pn) ratio. The DD F{P^a) generated by this 
diagram is just a constant |37|, which corresponds to a 
fiat TV = profile ho{/3, a) - 1/(1 - /3) and /(/3) -1-/3 
forward distribution. 



S{x-{kn)/{Pn)) 



ki = k + r/2 




r/2 



P + r/2 



P-r/2 



FIG. 4. Triangle diagram model for GPD. 

The calculation [35] of overlap integrals for light-front 
wave functions with a power-law behavior ilj{x^k±_) — 
l/{k'\_Y^'^ resulted in expressions equivalent to using 
DDs with N = Hi profile in Eq. ([29]) and forward dis- 
tributions behaving like (1 — The same profile 
arises [35 if one differentiates a scalar triangle diagram 
hz times with respect to masses (squared) of each active 
quark. 

The triangle diagrams, however, do not generate the 
Regge /(/3) — l/f3^ behavior for small /3. The latter 
may be obtained, in particular, by infinite summation of 
higher-order t-channel ladder diagrams (see, e.g., [38 ). A 
simpler way was proposed in Ref. [T5^ , where the specta- 
tor propagator was substituted by a parton-hadron scat- 
tering amplitude T(P, r, k) (see Figj5| written in the dis- 
persion relation representation. To avoid divergencies 
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ki 




r/2 



P + r/2 



P-r/2 



FIG. 5. Hadron-quark scattering amplitude. 



generated by the Regge behavior, the subtracted disper- 
sion relation 



functions that bring in an extra fall-oflF of the 

integrand at large transverse momenta k_\_. The authors 
of Ref. [15] intended to reflect this physics in their covari- 
ant model. To introduce form factors bringing in a faster 
fall-off of the /c-integrand with respect to quark virtual- 
ities kf and /c^, it was proposed to use higher powers of 
l/(m| — kf) instead of perturbative propagators, which 
may be achieved by differentiating the triangle diagram 
with respect to mf . 

The model of Ref. p!5] assumes spin-1/2 quarks. It 
was argued that the Dirac structure of the hadron-parton 
scattering amplitude in this case should be given by 
which provides EM gauge invar iance of the DVCS ampli- 
tude. Summarizing, the model scattering amplitude has 
the following structure 



T(P,r,fc) ^T{{P-kf)=To 



d(jp{(j) 



a-{P- fc)2 



(31) 



was used. The spectral function p(<j) here should be 
adjusted to produce a desired Regge- type behavior with 
respect to s = {P - k)'^ [39 . 

In the light-front formalism, the starting contribution 
corresponds to a triangle diagram in which the hadron- 
quark vertices are substituted by the light-front wave 



^ T{{P-kf 



(ml 



(32) 



To treat the two quarks on equal footing, it makes sense 
to take mi = m2, but we will keep them different for a 
while, to separate effects produced by nonzero A^i and 
N2. 

After ^ is contracted with the iji factor from the oper- 
ator vertex, one gets (/cn), and the model GPD that will 
be analyzed below is given by 



n^Ni+N2)\J (Pn) [m? - (A: + r)2]Wi+i[m2 - (A:-r)2]W2+i + '"^P^'^> \a - {P - 



a-{P-kf cr. 



(33) 



The overall factors were introduced here for future convenience. Using the a-representation 



I e^^^^-'^^^a^ da , -=/ e'^''' da^^ (34) 



(m2 - P)^+i 

for propagators and also for the l/cr subtraction term gives 

0^3 + {(^2 - ai)(rn)/2{Pn) 



POO nlOO 

— I da p{a) / ai^daia2'^da2das 
) ' Jo Jo 



S X 



{Ni^N2)\Jo Jo \ \ ai+a2 + a3 / (ai + 0^2 + as)^ 

(a2 - ai){rn)/{Pn) 



for the terms involved in the dispersion integral. The second delta-function corresponds to the 1/a subtraction term of 
the dispersion representation. It is accompanied by the l/(ai +0^2)^ factor because 1/a does not have /c- dependence. 
Introducing the skewness variable ^ = (rn)/2(Pn), changing = XiX and integrating over A we obtain 

x^^dxix^^dx2dx^5{l-xi-X2-x^) ( S {x - {x2 - xi)^) \ 
Taking equal masses mi = m2 = m, using xi -\- X2 = {I — xs) and introducing z through xi = (1 — xs)z results in 
"/o ^-P^-^l^^^^' [x,a + {1 - xs)mT^^^^^' [S{x - xs - (x, - x,)^) (i _ ,3)2 | - (37) 
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The To subtraction term gives the T)-term-type contribution 

< 1 ) (38) 



To ( x \ ( x\^W x'^^^ 



that vanishes outside the central region and, hence, is invisible in the forward limit. In what follows, we will concentrate 
on the terms generated by the dispersion integral, but one should remember that the I^o term can always be added 
to GPD i.e., in all formulas below one should be ready to change H{x,^) H{x,^) + Do{x/^). 



B. Forward case 



The case ^ = corresponds to the forward distribution 



H{x,^ = 0)=xj^ dapia) dx.dz _ ^^^^^^^^^^^ |<5(^ " ^s) - j^^, \ ■ (39) 

Taking x5{x) = for x 7^ gives 

"^'"''^ {N,+N2 + l)\ ^ ' Jo (a;(T + (l-a;)m2)^i+^2+i ■ 

Formally, we may also write 

'^(2^)^ dx^ZTTT^^^ (41) 



(40) 



Hix,,^0) ^m_,^^^ /,,3 fa) 

X X Jq X3{1-X3)^ 



The second term provides the subtraction regularizing the function f{x)/x at its singular point x = 0. 

C. DD description 



In the double distribution representation, we have x = /3 + a^. So, turning back to Eq. (36) and changing there 
1 — xi — X2 = /3, X2 — = a, we obtain that 

xi = ^{l-P-a) , X2 = + , (42) 

which gives for equal masses 

X r ^ f\o f^~\ (l-/?-a)^i(l-/? + a)^^ ^ ^, S{x-aO\ 



2A^i+A^2 + l 



Thus, a faster decrease of the /c-integrand with respect to the quark virtualities k\ or results in a suppression of 
the DD behavior by powers of (1 — /3 — a) or (1 — /3 + a) when a approaches the support boundary \a\ = 1 — j3. For 
equal Ni = N2 = we obtain 



22N+1 



Using Eq.(40), one can substitute the cr-integral through forward distribution to get 

X {2N+l)\ [\r^-l^ ^ [{^- -o?]"" f{P) \ ^ 5{x-ai,. 
«^M = ^,'-j^l dfi j da lf[six-P-aO- } . (45) 



This trick allows one to by-pass the question about the specific form of the spectral density p{cr). 
It is easy to notice that the factor 
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is precisely a normalized profile satisfying 

.1-/3 



hN{P,Oi)da = 1 . (47) 

Since xi -\- X2 < I for the Feynman parameters xi,X2, we have /3 > in the expressions above. The /3 < part of 
DD comes from the crossed diagram, in which the dispersion relation is written for T((P + /c)^). For the singlet case, 
the full DD should be symmetric with respect to interchange P —j3 (and also symmetric under a — Q^), 

which results in GPD i^(x,^) that is an odd function of x. For this reason, we will proceed with the /3 > case, 
keeping in mind to antisymmetrize the resulting H{x^ ^) at the very end. 
Thus, we can rewrite Eq.(45) as 



The first term 



^^^^^li^^ f\p r \a^hN{f3,a)6{x-P-a^) (49) 
X Jo J-i+(3 P 

coincides with the factorized DD Ansatz for H{x^ 0/^ which it is reconstructed from its forward limit f{x)/x. The 
relevant double distribution is given by f{P^a) = hN{P^ (^)f{P)/ f^- The second term may be rewritten as 



= -/ d(3 daS{x- 6{(3) 

Jo J-l-\-(3 Jo 



/ / . .,.^...^/ ^7(P^, (50) 



H{x,0 



and it provides a regularization of the /3-integral in Eq. 
(45). The total contribution is given by 

1 1^1-/3 

dp I daS {x — P — a^) 
J-i-\-(3 

xi^f{(3,a)-S{(3) '"'^7^^^! ■ (51) 

Thus, the model of Ref. [15], first, corresponds to the 
single-DD representation (22), and, second, it has the 



structure of the factorized DD Ansatz (27). Furthermore, 



due to the subtraction in the dispersion relation (31 ), one 
deals with the regularized double distribution 

r^iP, a) = fiP, a) - S{P) dj . (52) 

Returning back to Eq.(|48]) and calculating integral over 
a, we formally obtain 



H{x,O = ^0ix> 



(■Pi 

JI3i 



/(/?) 



h{/3, {x - /3)/0 



e{\x\ < 



JO 



d/3 



/3(l-/3)' 



hiP,x/0 , (53) 



where (see Fig. [6]) 



Pi — 7 > P2 



1-e 



(54) 




P 



FIG. 6. Support region and integration lines producing 
H{x,^) for X > <f and x < ^ from a double distribution that 
is nonzero for P > only. 



One should realize, however, that the second and the 
third integrals diverge and should be combined together 
to regularize their singularity at /3 = 0. This may be 
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achieved by rewriting the |a;| < ^ part in the form 



we obtain, for a; > ^: 



m 



/(/?) 



[h{p,{x-p)/0-hi/3,x/0] 
1 



3 X 



^2 



dp 



X- p 



explicitly showing the compensation of the 1/ j3 factor 

D. Results 

1. N = 1 profile 
For the model forward distribution 

= (1 - 

and the profile function 

3 (l-/3)2-a2 



hi{P,a) 



(58) 



4 (1 



Calculating i.e., the GPD at the border point 

X = one gets here the [(1 - - (1 - ^ P 

(55) 

factor from the profile function, and this factor changes 
the strength of singularity for /3 = 0. As a result, the 
integral over /3 converges as far as a < 1. This outcome 
is a consequence of using a profile function that linearly 
vanishes at the sides of the support rhombus. In its turn, 
the N = 1 profile is generated by the assumed l/(/^i/^2)^ 
dependence of the /c-integrand for large parton virtual- 
ities. If one takes the N = profile, the factor in the 
curly brackets should be substituted by 1/(1 — and 
the integral producing H{^^ ^) diverges. For small, but 
nonzero x — ^, one obtains the behavior proportional to 
l//3i {x — (,)~^. This result is similar to that obtained 
in Ref. [15 . However, since its authors explicitly declared 
that they are going to soften the hadron-quark vertices 
(56) by differentiating the diagram over the quark masses, one 
may wonder, how did it happen that they obtained a sin- 
gular result? 

The subtlety is that they took equal quark masses 
rui = 7X12 = m from the very beginning, and used dif- 
(^57) ferentiation with respect to this common m? . Here it 
should be noted that because 



^dm? ) (m^ — /c^)(m^ — /c^) (m^ — /c^)^(m^ — /c^) {m^ — k'lY^im? 



(59) 



the first and the third term on the r.h.s. are not softened with respect to one of the virtualities, i.e., one of the 
hadron-parton vertices remains pointlike. As we have seen above, imposing the l/(/c^)^^+^(/c2)^^^^ dependence on 
virtualities one would obtain the (1 — /3 — a)^^{l — /3 + a)^'^ /2^^~^^'^ factor, i.e., every differentiation with respect to 
ml gives (1 — /3 — a)/2, while every differentiation with respect to ml gives the (1 — /3 + a)/2 factor, both resulting in 
a nontrivial profile in the a-direction. On the other hand, each differentiation with respect to the common gives 
the {1 — P — a)/2 -\- {1 — P -\- a)/2 = {1 — P) factor that has no dependence on a. This kind of softening only increases 
the power of (1 — P), but DD remains flat in the a-direction. 



Note, that the use of l/(/cf )^+^(/c2)^^^-dependence in the model I^o-term contribution (38) results in the 
(1 — x^/^^)^ factor, which gives Do{l) = for A/" > case. This vanishing of Do{a) at the end-points a = ±1 
has the same nature as the vanishing of the DDs at the sides of the support rhombus: both result from a faster than 
perturbative decrease of the /c-integrand at large quark virtualities. 

Turning now to the |x| < ^ region, we use Eq. (55) to represent the relevant term for the A" = 1 profile as 



H{x,0\\x\<^ = 



3 X 

4 e 



(60) 



Note that as far as \x\ is strictly less than ^, the profile function does not vanish at the singularity point /3 = 0. 
The mechanism of softening singularity to 1/ p^ strength is now provided by the 1/cr subtraction term of the original 
dispersion relation. 

To get a model for singlet GPDs, one should take the antisymmetric combination 



Hs{x,0 = H{x,0-H{-x,0 . 



(61) 
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FIG. 7. Model singlet GPD Hs{x,^) with AT = 1 DD profile for a = 0.5 and ^ = 0.05,0.1,0.15,0.2,0.25. 

The resulting GPDs are shown in Fig. [2] For positive they are peaking at x = ^. The functions Hs{x,^) in this 
model are continuous at x = but the derivative dHs{x^£^)/dx is discontinuous at these points. 



2. N = 2 Profile 



Let us now take the N = 2 profile function 



and the same model forward distribution 



For X > ^ this gives 



15 X 
16^ 



15 [{1 - (3)^ - a^]^ 

16 (1-/3)5 

(1-/3)3 



/3^ 



./3)2 



X-/3 



(62) 



(63) 



(64) 



Evidently, the N = 2 profile gives a ~ /3^ suppression, and i^(^, ^) is finite as far as a < 2. 
Again, using Eq. (55), the |x| < ^ term can be represented in the form 



_i5x r^-i^i/^ dp 



+ (a; - Pf 
(1-/9)2 



15 X p d(3 



(/3-2) 



16 e 



^2(1-/3)2 



(65) 



explicitly showing the cancellation of the l//3^+^ singularity. 

The resulting GPDs are shown in Fig. [8) For positive x, they are peaking at points close to x = ^. In the model 
with N = 2 profile, both the functions Hs{x,(,) and their derivatives dHs{x^(^)/dx are continuous at x = 



E. D-term 



The 5 {13) subtraction term in the regularized DD (extended now onto the whole support rhombus) 

-l+|a| 



-l+|a| U 



1? 



(66) 
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FIG. 8. Model singlet GPD H{x, with N = 2 profile for a = 0.5 and ^ = 0.1, 0.3, 0.5, 0.7, 0.9. 



softens the singularity of a) for /3 = 0, but it does not convert a) into a "plus distribution" a)]+ whose 
integral over /3 vanishes. Thus, contains a nonzero D-term contribution 



r^{p, a)d/3 = a dp ^ h{f3, a) { 1 

-l+|a| i-l+|a| 



l-|a| B-l 

dp m HP, a) p 



(1-/3)2 ■ 
(67) 



Taking the same model forward distribution = (1 — j3Y / and = 1 profile function gives 



^ ' 2 



(1 



{(3-2). 



(68) 



A similar expression for the i^-term is obtained in the N = 2 profile model: 



^ ' S Jo 13- 



a 



(1-/5)'J 



2). 



(69) 



As one can see in Fig. [gj the two curves are rather close 
to each other. 



D{a) 




FIG. 9. The D-terms in AT = 1 and = 2 profile models for 
a = 0.5. 



The comparison of the total GPD and its 

D-term part is shown in Fig. [lO] The difference between 
GPD and D-term D{x/(^) corresponds to the 

term obtained from the "plus" part + 

of DD. The shape of the difference for ^ = 0.5 is shown 




FIG. 10. GPD and 7:>-term D{x/^) for ^ = 0.5 and 

positive X. 



in Fig. [TT] Note that, despite the fact that the forward 
distribution in this model is positive, there is a region, 
where the contribution to H{x, ^) coming from + 
is negative. This is due to the S{/3) subtraction term con- 
tained in + . Also shown is the ratio H_^{x,^)/x. 
Looking at the figure, one may suspect that the x-integral 
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FIG. 11. Difference between GPD H{x,^) and D-ieim. 
D{x/^) in the case of the = 1 profile for ^ = 0.5 and 
positive X. The same function divided by x is also shown. 



of vanishes. In the next section, we show that 

this, indeed, is the case. 



can be analytically continued to the point j = —1. The 
result of such a procedure is equivalent to analytic reg- 
ularization of the x-integral. However, the required an- 
alyticity properties of ^{j) may be violated by singular 
or "invisible" terms (cf. [20 ) in the integrand of Eq.(72) 
(e.g., xS{x) gives a non-analytic Sj^-i contribution into 
^(j)). In the model construction described above, sin- 
gular terms explicitly appear as a result of subtractions 
in the dispersion relation, so our intention is to develop 
a less restrictive approach to this problem. 



Below, we give a derivation of the sum rule ( 71 ) based 
on separation (25) of the DDs into the "plus" part and 
the D-teim. No assumptions about smoothness will be 
made. In fact, the key element of the derivation is that 
H{x^x)/x should be treated as a (mathematical) distri- 
bution at the point x = rather than a function. The 
same applies to H{x^O)/x. 



IV. GPD SUM RULES 
A. Formulation 

The D-teim determines the subtraction constant in the 
dispersion relation for the DVCS amplitude p^H22] . In 
particular, it was shown [19 that the original expres- 
sion for the real part of the DVCS amplitude involving 
H{x, ^), and the dispersion integral involving H{x^ x) dif- 
fer by a constant A given by the integral of the D-teim 
function D{a): 



J-i x-^ J_i 1 - 



a 



Here, P denotes the principal value prescription. In 
Ref. [19], this relation was derived using polynomiality 
properties of GPDs. It was also pointed out there that it 
can be obtained by incorporating representation of GPDs 
in the two-DD formalism (which is basically again the use 
of the polynomiality). 

Taking ^ = 0, one formally arrives at the sum rule 



/: 



^ H{x,x) - H{x,0) 



dx 



-fj 



' D{a) 



da . 



a 



(71) 



Since both H{x^O)/x and H{x^x)/x are even functions 
of X, their singularities for x = cannot be regular- 
ized by the principle value prescription. Moreover, there 
are no indications that singularities of these two func- 
tions may cancel each other. On the contrary, as em- 
phasized in Ref. [40 , there are arguments that the ratio 
H{x^x)/H{x^O) does not tend to 1 for small x. 

The solution given in Refs. [20l ISH |4T] is based on the 
analytic regularization of the x-integral. Namely, it is 
assumed that the positive Mellin moments (or conformal 
moments, see, e.g., [42^ ) 



m 



/: 



x^[H{x,x)-H{x,0)] dx 



(72) 



B. Ingredients 

To begin with, we remind the basic formulas: the ex- 
pression 



H{x,0 



/(/3, a) S{x-P- ^a) dp da (73) 



X JQ 

producing GPDs from DDs and the decomposition of DD 

Dia) 



/(/?,a) = [/(/?, a)]++5(/?): 



into the "plus" part given by 



[/(/?, a)] + = /(/?, a) -5(/?) f /(7,a)rf7 

J-l+\a\ 



(74) 



(75) 



and the D-teim part S{P)D{a)/a. 

Correspondingly, we split GPD into the part coming 
from the "plus" part of DD 



H+ix,0 _ 



S{x — P — ^a) — 5{x — ^a) 



and that generated by the D-ieiia 
Hd{x,0 _ D{a) 



S{x — ^a) da 



d/3 da 
(76) 

(77) 



X J_i a 

The latter integral gives an explicit expression 

Hd{x, = sign(0 0{\x\ < Id) D{x/0 , (78) 

but, as we will see, it is instructive to use the integral 
representation as well. Another important relation 



--S{x) 



f' D{a 

7-1 



da 



(79) 
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is obtained by taking ^ = 0. 

Note, that if we take x = ^, Eq. (78) gives 



Hd{x, x) = sign(x) D(l) 



(80) 



If D(l) ^ 0, then the integral for A in ([70]) di- 
verges. However, as argued in the previous section, 
1^(1) = for models with faster-than-perturbative de- 
crease of the hadron-parton amplitude at large quark 
virtualities. Thus, we assume that D{1) = 0, and fur- 
thermore that the integral of D{a)/{1 — a) converges. 
Then Eq. ^Tj) gives 



--6{x) 



a{l — a) 



da . 



(81) 



The important point is that if we would use this formula 
to write an expression for Hd{x^x) itself, we would get 
x5{x) on the r.h.s., which should be treated as zero for 
integration with functions finite for x = 0, since the co- 
efficient given by the a-integral is also finite. Thus, the 
scenario with D(\) = is self-consistent. 

Note that both Hd{x^^) /x and H]j{x^x)/x are pro- 
portional to S{x), with the coefficients given by integrals 
of D{a). This means that, unlike the functions H{x^O) 
and which, for x 7^ 0, are insensitive to changes 

of D{a) in the S{P)D{a)/a term, the (mathematical) dis- 
tributions H{x^O)/x and H{x^x)/x contain information 
about such a D-teim. 

Our next step is to study contributions from different 
parts of the GPDs involved in the sum rule (71). 



'Secondary" sum rule 



1. "Plus" part 



a. Forward function One can easily see from 
Eq. ([76|) that 



dx = 



(82) 



for any ^, including ^ = 0. Since the integrand is an even 
function of x, the vanishing of this integral means that 
we also have 



dx = . 



(83) 



Thus, i^+(x,^) should be negative in some part of the 
central region, and this negative contribution should 
exactly compensate the contribution from the regions, 
where i^+(x,^) is positive. In other words, on the (0, 1) 
interval, i^+(x,^)/x has the same property as a "plus 
distribution" with respect to x. Note, that this does not 
mean that H_^{x^^)/x necessarily contains singular func- 
tions like 5{x). For finite ^, the function H-^{x^^)/x is 



pretty regular for all x values (see Fig 12). The negative 



6{x) function appears only in the ^=0 limit, i.e. 



H+{x,0) fix) 



5{x) 



fix) 



(84) 



(Here, it was taken into account that i^+(x,0) coincides 
with the forward distribution f{x) for x 0). 




X 



0.4 0.6 0.8 1.0 



FIG. 12. Function H^{x,^)/x in the model with the = 1 
profile for ^ = 0.2, 0.3, 0.5 and positive x. 

b. Border function For the integral involving the 
border function, we get 



dx 



dx / dp da /(/3, a) 

1 



X [6 (x(l - a) 



8{x{l-a))\ . (85) 



Noting that equation a = 1 is satisfied in one point on 
the support region Vt only, namely, in the upper corner 
of the rhombus, we may treat 5{x{l — a)) as 5{x) /(I — a) 
to get 



r'^±^cix=f\x [ 

J-l X J_i 



<5 X 



1 — a 



5ix) 



■ (86) 



Since |/3/(l — a)| < 1, the first (5-function always works. 
As a result, the integrals coming from the two delta- 
functions cancel each other, and we have 



' H+ix,x) 



dx = ^ 



(87) 



just like for i7+(x,^)/x. Unlike i^+(x,^), however, the 
combination H_^{x^x)/x explicitly contains the S{x) sub- 
traction term, i.e. it is a genuine "plus distribution" with 
respect to x: 



Hj^{x^x) H{x^x) 



Six) 



dy 



H{x^ x) 
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Summarizing, the "plus" parts of both functions en- 



tering into the sum rule ( 71 ) separately produce vanish- 
ing contributions into the x-integral. Furthermore, these 
zero contributions are due to the fact that i^+(x,0)/x 
and Hj^{x^x) /x are "plus distributions", which results 
in zero integrals irrespectively of the form of the forward 
distribution fix) and the border function H{x^x). 



2. D-part 



Let us now turn to the D-paits. First, we have 



dx 



' D{a) 



da 



(89) 



for any fixed ^, including ^ = 0. This result may be 
obtained by integrating over x the S{x — ^a) factor in the 
integral representation (77). For non-vanising ^, one can 
also use Eq. ( jTS] ) in the x-integral and then change the 
integration variable through x = a£^. 

For the integral involving the border function, we use 
Eq. (81), which gives 



/: 



^ Hd{x,x) 



dx= 



D{a) 
J_i a{l - a) 



da . 



As a result. 



7-1 X J_i X J_i 1 



Dia) 



(90) 



da . 
(91) 



Combining this outcome with zero contributions from the 
"plus" parts, one obtains the sum rule (71). 

Thus, our construction confirms the sum rule. But 
our derivation shows that the "plus" parts of both terms 
simply do not contribute to the sum rule whatever the 
shapes of f{x) and H{x^x) are. Only the I)-parts con- 
tribute, so there is no surprise that the net result can be 
expressed in terms of D{a). 

An essential point is that both Hd{x,0)/x and 
Hjjix^x) / X are proportional to the (^(x)-function, with 
the coefficients given by integrals of the D-ieiia function 
D{a). In this sense, H{x,0)/x and H{x^x)/x "know" 
about the D-teim. 

A simple consequence is that all x^ moments of 
Hd{x^ 0) and H]j{x^ x) with j > vanish and one cannot 
get the I)-part of the sum rule fn\ by an analytic con- 
tinuation of the x^ moments of H]j{x^O) and Hd{x^x) 
to j = —1, i.e., using the procedure of Refs.[20l [SU |4T] . 
In fact, x^ moments of Hd{x^O) and Hd{x^x) are pro- 
portional to the Kronecker delta function ^7,-1. 



3. Formal derivation and need for renormalization 

Since H{x^ 0)/x is given by integrating the DD a) 
over a along vertical lines /3 = x, a subsequent integra- 
tion over all x gives DD f{P^a) integrated over the whole 



rhombus: 

H{x,0) 



f ^^^^dx= j dx j dpdaf{f3,a)S{ 



x-p) 



I f{l3,a) dl3da= f 
Jn J-i ct 



da . 
(92) 



On the last step, we used that the /3-integral of a) 
formally gives D{a)/a. However, if f{P^a) ~ 
being even in /3, one needs a regularization for the 
/3-integral. The "DD+D" separation ( [73| ), as we have 
seen, provides such a regularization. It works like a renor- 
malization: the divergent integral formally giving the D- 
term is subtracted from the "bare" DD, and substituted 
by a finite "observable" function D{a)/a. 

In a similar way, we can treat the second integral: 



^ H{x,x) 



dx= 



/dx / dp daf{l3^a)5{x — P — xa) 
-1 Jn 



f 11 

Jn 1 



dp da 



J_i a{l - a) 



da . 
(93) 



Again, the last step requires a subtraction of the infinite 
part of the /3-integral. 

The advantage of using the "DD+D" separation as a 
renormalization prescription is that it is applied directly 
to the DD. Hence, it is universal, and may be used for 
other integrals involving 



4. Comparison of the "plus" prescription and analytic 
regularization 

Another possibility to renormalize the /3-integral for 
a singular DD is to use the analytic regularization em- 
ployed in Refs.[20l [SH |4T]. It works as follows. If we 
need to integrate a function like X{x)/x^~^^ with X{x) be- 
ing finite and nonzero for x = 0, we subtract from X{x) as 
many terms of its Taylor expansion as needed to remove 
the divergence 



y x{x) ^ py ^ A(x)-A(o)-xy(o)-... 

dx = dx 



(0) 



v.a+1 



r dx fy dx 

J{0) ^ J{0) 



(94) 



and then treat the compensating integrals of x'^/x^'^^ as 
convergent, substituting them by y^~^/{n — a). 

So, let us consider again a DD which is nonzero for 
positive P only and has the form 

f(P,a) = ^^^0{/3+\a\<l)0{/3>O) 



16 



with a < 1. Then the analytic regularization of its inte- 
gral with some reference function $(/?) is defined by 



^(/?)A(/3,a) 



d/3 



(0) ^' 

/o /^"+^ 
which may be rewritten as 

.l-\a\ 



dp- 



(95) 

<^(Q)A(Q,a) 

a(l - \a\Y ' 



/ 

J(0) 



/o 
•$(0) 



Af/3 a) 



A(/3,a)-A(0,a) 



(96) 



a+1 



A(0,a) 
a(l - \a\Y 



Now, the first contribution on the r.h.s. is generated by 
the "plus" part of the DD, while the second one comes 
from a D-term. After adding the /3 < part of the 
DD, the D-ieiia V{a)/a corresponding to the analytic 
regularization is given by 



V{a) 



a 



A(;9,a)-A(0,a) 



/3' 



a+1 



dp 



A(0,a) 
a(l - lal)'^ 



(97) 



Thus, the analytic regularization prescription unambigu- 
ously fixes the D term, and in this sense it may be called 
the "analytic renormalization" . 

In the model considered in the previous section, we also 
obtained a concrete result for the D-ieiia. But the spe- 
cific i^-term contribution we obtained there came only 
from the cr-integral part of the dispersion relation for 
the hadron-parton scattering amplutude subtracted at 
{P — =0. As we pointed out, one should be always 
ready to add to it the I^o term coming from the Tq con- 
stant in the dispersion relation (31 ). In principle, we had 
no reasons to require that Tq = 0. In this sense, the 
D-teim in that model is not fixed. 

On the other hand, the statement, that moments 
of iJ(x,^) are analytic functions of j, does not explic- 
itly mention fixing any subtraction constants: it sounds 
like a general principle, and may create an impression 
that there are no ambiguities in the subtraction of the 
/3 = singularity. However, the analyticity assumption 
was not shown so far to be a consequence of general prin- 
ciples of quantum field theory. Moreover, as mentioned 
in Ref. [43] , it is not satisfied in the nonlocal chiral soliton 
model. Still, one may hope that it is valid in QCD. 

To see if the Tq = model of the previous section 
agrees with the analyticity assumption, we should just 
check whether its D-teim is different from that obtained 
via analytic renormalization. In particular, for the = 1 
model, we have 



a 



(98) 




FIG. 13. The D-terms in the model with A = 1 profile and 
a = 0.5: T>(a) was obtained using analytic regularization, and 
D{a) was obtained for To = in the model of the previous 
section. 



and, hence. 



a 



(1-H)2- 



,(l-|a|)i-« 



1 



o(l - \a\)°- 
(99) 



In Fig 13, we compare this result (for a = 0.5) with the 
result obtained by single subtraction in the dispersion 
relation dsTl) with Tq = 0. 



Our main point is that representing H{x^£^) as the sum 
JT+(x,^) + Hd{x^(_) one can derive the GPD sum rule 
( [7l| without using the analyticity assumption. But since 
our derivation, so to say, works for any D-term, it also 
works for the D-teim following from the analyticity as- 
sumption. 



D. Generic sum rule 

Finally, let us apply the "DD+D" separation to the 
generic relation ( 70 ) . 



1. "Plus" part 



Representing 



1 



X — ^ X {x — ^)x 



(100) 



and using Eqs. ( |83[ ), ( |87[ ), we have 

dx 

dx = H I F - 

-1 x-^ 



1 x-^ 

5{x{l-a) - p) -5{x{l-a)) 



I /(/3, a) dp da 
(101) 



Jn 



) dp da 



p-^{l-a) {I -a) 
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and 



P f ^^±^dx = P [\—^ [ f{/3,a)dpda 
J-1 ^ — ^ J-1 X — ^ Jn 



6{x — P — ^a) — 5{x — ^a) 



(102) 



) dp da 



p-^{l-a) {I -a) 



Thus, seemingly different delta-functions have converted 
— ^) into identical expressions (cf. Ref. [21 , where a 
similar result was obtained for the Fjj part of the two-DD 
representation). As a result, 



J-1 x-^ J_i 



■dx-P I -^^^^y dx = . (103) 

In this case, we deal with the situation when the differ- 
ence of two integrals vanishes, but each integral does not 
necessarily vanish. 



2. ''Depart 

For the integral involving the border function, we have 

H]j{x^x) 



dx 



dx — 
1 X- 



S{x) 



1 X 

a{l — a) 



x-i 
da = . (104) 



In simple words, the starting integrand in (104) vanishes 



for X 7^ since then Hd{x,x) = 0, while for x = it 
is given by the xS{x) distribution which produces zero 
after integration with a function that is finite for x = 0, 
which is the case if ^ 7^ 0. Comparing this result with 
Eq. ( 90 ) , we see that the nonzero value given by the latter 
cannot be obtained by taking ^ = in the final result of 
Eq. (IToa above. 



The second piece is given by 



J-1 x-^ J_i X X 

^ xdx D(a) ^, ^ ^ -, 
P / / S(x - ^a) da 

J-ix-i y_i a 

ia D{a) _ D{a) 



dx 



da 



da 



(105) 



Again, the result above may be obtained by simply us- 
ing Hd{x,(_) = sign(^) 6>(|x| < \(,\)D{x/(_) and rescaling 
X = a(_. Also, though the final result of Eq. ( 105 ) does 
not depend on ^, it does not coincide with the result of 
the counterpart relation (89). 



However, for the difference of the two integrals we ob- 
tain 



Hpjx.x) 
-1 x-^ 



dx 



' D{a) 



da 
a 

(106) 



the same result as in Eq. (91). Combining the results for 
the "plus" and I)-parts gives Eq. (70). 



3. Some conclusions 

Thus, our calculation confirms the generic GPD sum 
rule ([70| derived in Refs. [19, 21j. We were also able to 
derive the ^ = sum rule (71) suggested in Ref. [19 . It 



should be emphasized that the integrals present in the 
generic sum rule have a singularity for x = ^, which 
is inside the region of integration, so the integrals may 
be taken using the principal value prescription. Since 
H{x^O)/x and H{x^x)/x are even functions of x, the 
^ = sum rule may be written through an integral from 
to 1, and its 1/x singularity is at the end-point of the 
integration region, which means that the P-prescription 
cannot regulate it. Just because of this fact alone, the 
sum rule (71) cannot be a straightforward consequence 



of the generic sum rule (70). 



In our derivation, we managed to obtain finite expres- 
sions for each term involved. In particular, we established 
that though Hd{x^x) and Hb{x^£) contributions to the 
generic sum rule (70) are ^-independent, they do not co- 



incide with their counterparts from the secondary sum 
rule (71), i.e., the latter cannot be obtained by formally 
continuing to ^ = the ^-independent results for each 
term of the generic GPD sum rule. 

In our derivation, we did not make an assumption 
about analyticity of the Mellin moments of GPDs. We 
have obtained GPD sum rules as a consequence of the 
polynomiality of GPDs that follows from Lorentz invari- 
ance and is encoded in the DD representation. The an- 
alyticity is a much stronger restriction. One may try to 
find out whether it can be tested experimentally and it 
is also worth trying to prove it in QCD. 



V. SUMMARY 

In this paper, we discussed some basic aspects of build- 
ing models for GPDs using the factorized DD Ansatz 
(FDDA) within the "single-DD" formulation. The main 
difficulty in the implementation of such a construction is 
the necessity to deal with projection onto a more singu- 
lar function f{/3)//3 (rather than just onto /(/3)) in the 
forward limit. This leads to two problems. First, one 
encounters non-integrable singularities for /3 = in the 
integrals producing GPDs in the central region \x\ < |^|. 
The difficulty is exaggerated by necessity to consider for- 
ward distributions f{f3) that have a singular Regge 
behavior at small p. Second, if there are no factors sup- 
pressing the P ^ region for the integration line cor- 
responding to X = ^, the combined l//3^+^ singularity 
leads to a singular {x — ^)~^ behavior for GPDs in the 
outer region x > ^ near the border point x = ^. Such a 
behavior was found in the model of Ref. 15 . 
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In our analysis, we found that this model gives the 
single-DD-type representation for the model GPD, and 
thus above reasoning is applicable to it. But we argued, 
that a proper softening of the hadron-quark vertices pro- 
duces a profile function /i7v(/3,a) that results, for x = ^, 
in the suppression factor securing a finite value 

of the GPD at the border point. 

However, the profile factor has no impact on the com- 
bined l//3^+^ singularity on the /3 = line inside the sup- 
port rhombus, which one faces when calculating GPDs 
in the \x\ < |^| region. The advantage of the model of 
Ref. [15 is that it implants the Regge behavior through a 
subtracted dispersion relation for the hadron-quark scat- 
tering amplitude. We found that the subtraction pro- 
vides the regularization necessary for the calculation of 
GPDs in the central region, and illustrated the behavior 
of resulting GPDs in models with = 1 and N = 2 
profiles. 

We also observed that this model produces a D-teim 
contribution, despite the fact that it uses only the for- 
ward distribution as an input. This D-teim contribution 
appears because the subtraction generated by the disper- 
sion relation differs from the subtraction that converts 
the original DD into a "plus" distribution + . 

The latter, by definition, cannot generate a D-teim. We 
have shown that the GPD generated by the 

part of the original DD (i.e., GPD H{x,^) 
with the D-teim contribution D{x/(^) subtracted) has a 
remarkable property that the integral of i^+(x, ^)/x over 
positive values < x < 1 vanishes. As a result, iJ+(x,^) 
must be negative in some part of the central region, a 
feature that is absent in previous FDDA models based 
on two-DD formulation. 

Within the single-DD formalism, it is very natural to 
separate the relevant DD /(/3,a) into the "plus" part 
[/(/3,a)]+ and the D-ieiia. We demonstrated that this 
separation can be used to rederive the GPD sum rule 
related to the dispersion relation for the real part of 
the DVCS amplitude, and we also gave a derivation of 
another sum rule proposed as the ^ ^ limit of that 
generic sum rule. Our derivation shows that this "sec- 
ondary" sum rule is not a straightforward consequence 



of the generic one. In particular, the principal value pre- 
scription used in the generic sum rule needs to be substi- 
tuted by another prescription, like the "plus" prescrip- 
tion. The "plus" prescription, in fact, is automatically 
generated by the separation of DDs into the "plus" part 
and the D-ierm. We also demonstrated that the contri- 
butions into the two sum rules generated by the same 
functions are not in a one-to-one correspondence. 

Summarizing, using (intentionally) simplified models, 
we developed the basic tools that can be used in build- 
ing realistic GPD models based on the factorized DD 
Ansatz within the single-DD formalism. Future develop- 
ments in this direction should include the extension of 
the presented methods onto the cases with a > 1 Regge 
behavior, which would require an extra subtraction in 
the dispersion relation, and building models for nucleons 
and other targets with a non-zero spin. 



ACKNOWLEDGMENTS 

I would like to express my deep gratitude to A. P. 
Szczepaniak for numerous and intensive communica- 
tions about Ref. [15] that initiated this work. I thank 
D. Miiller, M.V. Polyakov, K.M. Semenov-Tian-Shansky 
and O.V. Teryaev for stimulating discussions of GPD 
sum rules. I am grateful to I.V. Anikin, I.I. Balitsky, 
A.V. Belitsky, S.J. Brodsky, M. Burkardt, M. Diehl, 
M. Guidal, V. Guzey, C. E. Hyde, C.-R. Ji, X. D. Ji, 
P. Kroh, S. Liuti, J.A.Miller, I.V. Musatov, M.V. 
Polyakov, A.Schafer, M.A. Strikman, L.Szymanowski, 
A.W.Thomas, B. C. Tiburzi, M. Vanderhaeghen and 
C. Weiss for many inspiring discussions and communi- 
cations that we had over the years, and which eventually 
infiuenced this work. 

This paper is authored by Jefferson Science Asso- 
ciates, LLC under U.S. DOE Contract No. DE-AC05- 
06OR23177. The U.S. Government retains a non- 
exclusive, paid-up, irrevocable, world-wide license to 
publish or reproduce this manuscript for U.S. Govern- 
ment purposes. 



[1] D. Mueller, D. Robaschik, B. Geyer, F. M. Dittes, and 
J. Horejsi, Fortschr. Phys., 42, 101 (1994), ,arXiv:hep^ 

l^ph/9812448 

[2] X.-D. Ji, Phys. Rev. Lett., 78, 610 (1997)| |arXiv:hep-| 

[ ph/9603249l 

[3] A. V. Radyushki n, |Phys. Lett., B380, 417 (1996)1 
arXiv:hep-ph /96043T7I 



[4] J. C. Collins, L. Frankf urt, and M. Strikma n, [Phys. Rev., 
[ D56, 298 2 (1997), ,arX iv:hep-ph/9611433| 

[5] X.-D. Ji, I J. Phys., G24, 1181 (199^)] [ arXivihep- 
[ ph/9807358 

[6] A. D. Martin and M. G. Ryskin, ^Ph ys. Rev., D57, 6692 
[ (1998), arXiv:hep-ph/9711371 

[7] B. Pire, J. Soffer, and O. Teryaev, |Eur. Phys. J., C8, 



103 (1999), arXiv:hep-ph/9804287^ 

[8J A. V. Radyushkin, Phys. Rev., D59, 014030 (1999)| 



[9] 
[10] 

[11] 
[12] 

[13] 



[a?Xiv:hep-ph/9805342 
A. V. Radyushkin, | 



|Phys. Lett., B385, 333 (1996)| 



|arXiv:hep-ph/9605431 
"Dual parameterization" [401 1411 1431446] is another way 
to impose the polyno miality condition onto model GPDs. 
A. V. Radyushkin, jPhys. Lett., B449, 81 (1999) 
arXiv:hep-ph/9810466 

M. V. Polyakov and C. Weiss, |Phys. Rev., D60, 114017J 

(1999), arXiv:hep-ph/9902451 

K. Goeke, M. V. Polyakov, and M. Vanderhaeghen, 
JProg. Part. Nucl. Phys., 47, 401 (2001)[ |arXiv:hep-l 
ph/0106012j 



19 



A. V. Belitsky, D. Mueller, A. Kirchner, and A. Schafer, 
Phys. Rev., D 64, 116002 (2001) , arXiv:liep-pli/QQ11314 . 
A. P. Szczepaniak, J. T. Londergan, and F. J. 
Llanes- Estrada, Acta Phys. Polon., B40, 2193 (2009), 
arXiv:0707.1239 [hep-ph] 

P. V. Landshoff, J. C. Polkinghorne, and R. D. Short, 
Nucl. Phys., B28, 225 (1971) 

K. Kumericki, D. Mueller, and K. Passek-Kumericki, 
(2007), arXiv:0710.5649 [hep-ph] 
O. V. Teryaev, (2005), arXiv:hep-ph /05 100311 



I. V. Anikin and 0. V. Teryaev, Phys. Rev., D76, 056007 


(2007), arXiv:0704.2185 [hep-ph] 




K. Kumericki, D. Mueller, anc 


K. Passek-Kumericki, 


Nucl. Phys., B794, 244 (2008) 


arXiv:hep-ph/0703179 


M. Diehl and D. Y. Ivanov, |] 


Eur. Phys. J., C52, 9191 


(2007), arXiv:0707.0351 [hep-ph^ 





O. Teryaev, PoS, DIS2010, 250 (2 010). 
X.-D. Ji, JPhys. Rev., D55, 7114 (1997)| |arXiv:hep- 
ph/960938ll 



[30] 



A. V. Radyushkin, (2000), arXiv:hep-ph/010122^ 

M. Diehl, Phys. Rept., 388, 41 (2003)^ ,arXiv:hep- 

ph/030738l] ~_ 

A. V. Belitsky and A. V. Radyushkin, |Phys. Rept., 418 

1 (2005), arXiv:hep-ph/0504030 

S. Boffi and B. Pasquini, Riv. N uovo Cim., 30, 387| 

(2007) , arXiv:0711.2625 [hep-ph] 

A. V. R adyushkin, Phys. Rev. , D58, 114008 (19"98)j 
■arXiv:hep-ph/9803316 

DDs and GPDs depend on the momentum transfer 
t = r^, but this dependence is not important for our pur- 
poses. So, in what follows, we consider the formal t = 
limit. 

Here and below we consider, for simplicity, spin-0 
hadrons. 



[35 



[36 
[37 
[38; 
^39 

[40 



I. I. Balitsky and V. M. Braun, [Nucl. Phys., B311, 541 

(1989) 

O. V. Teryae v, [Phys. Lett., B510, 125 (2001)||arXiv:hep-| 
ph/010 2303' 

B. C. Tiburzi, iPhys. Rev., D 70, 057504 (2004)] 
arXiv:hep-ph/0405211 

K. Kumericki, D. Mueller, and K. Passek-Kumericki, 
|Eur. Phys. J., C58, 193 (2008)[ |arXiv:0805.0152 [hep-] 

A. Mukherj ee, I. V. Musatov, H. C. Pauh, and A. V. 
Radyushkin, [Phys. Rev., D67, 073014 (2003)| |arXiv:hep- 

ph/0205315 

D. S. Hwang and D. Mueller, Ph ys. Lett., B660, 350] 

(2008)^, arXiv:0710.15 67 [hep-ph] 

A. V. Radyushkin, JPhys. Rev., D56, 5524 (1997) 



arXiv:hep-ph/9704207| 

A. Efremov and A. Radyu shkin, [Mod. Phys. Lett., A24,J 
2803 (2009), arXiv:0911.1195 [hep-ph] 
To get the s"' Regge behavior with 1 < a < 2, one should 
use a doubly subtracted dispersion relation, but in this 
paper we will follow the original construction of Ref . [15] , 
leaving the generalization for a > 1 to a future work. 
M. V. Polyakov, Phys. Lett., B659, 542 (2008) | 
|arXiv:0707.2509 [hep-ph] 



M. V. Polyakov and K. M. Semenov-Tian-Shansky, Eur.| 
Phy s. J., A40, 181 (2009), arXiv:0811.2901 [hep-ph] 
K. Kumericki and D. Mueller, [Nucl. Phy s., B84l7T] 



303 



_(2010)",'arXiv:0904.0458 [hep-ph] 

J<:. M. Semenov-Tian-Shansky, Eur. Phys. J., A36 
(2008) , arXiv:0803.2218 [hep-ph] " 
M. V. Polyakov and A. G. Shuvaev, (2002), [arXivihep- 
_ph/0207153 

M. V. Polyakov, (2007), arXi v:071 1.1820 [hep-ph]] 



K. M. Semenov-Tian-Shansk y, |Eur. Phys. J., A45, 217 
(2010)narXiv:1001.2711 [hep-ghj] 



